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ABSTRACT 


The main objective of this study is to find a wavelet-based, feature extracting 
algorithm for push-to-talk transmitter identification. A distance-measure algorithm 
is introduced to classify signals belonging to one of four transmitters. The signals 
are first preprocessed to put them into a form suitable for wavelet analysis. The 
preprocessing scheme includes taking envelopes and differentials. Median filtering 
is also applied to the envelopes and the differentials to denoise the data. The pre- 
processed data takes on a pulse-like shape, which is suitable for wavelet processing. 
The distance algorithm is applied to the outputs (scales) of the wavelet transform. 
The distance algorithm uses local extrema of the wavelet coefficients, and computes 
the distance between the local extrema of a template and the processed signals. A 
small distance implies high similarity. A cone from each transmitter is selected as 
a template. The distance measure is computed between any signal of interest and 
the reference templates. The signals are identified to belong to one of four transmit- 
ters according to the distance measure. A small distance measure indicates that the 
signal belongs to the transmitter from which the template originated. The distance 
algorithm can classify correctly the four different signal sets provided for the research 


and, even at lower signal-to-noise levels, good identification is achieved. 
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I. INTRODUCTION 


Signals produced by different transmitters have different characteristics. The 
transient response, obtained when the transmitter is turned on or off, is basically a 
build-up from off-state to on-state, or vice versa. These responses are unique to the 
specific transmitter and differ even among units of the same make and model. This 
observation applies only to the turn-on response. Therefore, a technique to extract 
the uniqueness of the turn-on transients will help to identify the transmitter. 

Time-frequency analysis of stationary signals is a well-known subject. The 
Fourier Transform method is quite efficient for this kind of analysis. However, when 
the signal of interest is non-stationary, the Fourier Transform method is not appro- 
priate, because it uses as a basis function the complex exponential that expands over 
infinite time. That is, it lacks the time-evolution of frequencies that is very important 
for the analysis of transients of short duration. 

A sliding time-window was introduced by Gabor (1946) to gain time informa- 
tion from the Fourier Transform method. This modified Fourier Transform method 
is called the Short-Time Fourier Transform (STFT), and uses a modulated window 
as a basis function. The STFT looks through the window, and processes just a small 
portion of the signal and, by sliding the window, processes another small portion of 
the signal successively. Thus, once a window has been chosen, time-frequency reso- 
lution is fixed. The problem with this method is that the signal must be stationary 
within these small portions, or else the method will have limitations in reflecting the. 
time-evolution of frequencies within the window. 

Another technique for transient Sena analysis is the Wavelet Transform (WT), 
which seems more efficient in terms of time evolution of frequencies than the STFT. 


Basis functions of the WT, unlike the complex exponential of the Fourier Transform, 





do not have infinite duration. They are nonzero for only a short duration. Accord- 
ingly, this “compact support” makes the WT localized, not only in frequency but 
in time as well. Moreover, wavelets provide the flexibility to choose the particu- 
lar wavelet function that is appropriate for a specific application. This is possible 
since there are a large number of compactly supported wavelets that can be used as 
orthogonal basis functions. 

The purpose of this thesis is to investigate the use of the Wavelet Analysis 
for the classification of transient signals of push-to-talk transmitters, and to find an 
algorithm to identify transmitters. Chapter II gives a brief introduction to wavelet 
analysis and median filtering. Chapter III explains the algorithms to identify the 
push-to-talk transmitters, and Chapter IV shows the results and presents the recom- 


mendations for future studies. 








Il. SIGNAL DECOMPOSITION AND MEDIAN 
FILTERING 


A. CONCEPT OF A BASIS 


A linear independent set S of vectors {v,} forms a basis of the space V if 


every element x € V can be represented in the form of a series 


x = S- CeVrE , (21) 
k 


where c;’s are the coefficients. 
A set S of vectors are said to be linearly independent if and only if there are 


numbers ¢),C9,..-Cn,---, all of which are zeros, such that 


CyVy t+ CoVoa t+... + CnVn t... = 0. 


B. INNER PRODUCT SPACES 


An inner product space is a vector space V with an inner product defined 
on it. Ifx = (21,29,...,2n) and y = (y1, Yo,---, Yn) are vectors in V, then the inner 


product of x and y are defined by 


xy) =x?.y* =a yf + a.yf+..+2n.y% (272) 


Thus, an inner product of x and y is a mapping from x.y, cartesian product space, 
into the scalar product field. 
The inner product has the following properties: 


Let x, y, z be vectors in V, and c any real-valued scalar, then 
(xX + y),Z) = (x2)+( yz). 
x,x) >0O with equality if and only ifx = 0. 


Functions like f (xz), on the interval 0 < x < X, can be thought of as a vector 


with a continuum of components along the whole interval. In this case, the summation 


3 


can be replaced by integration (Burrus and Gopinath, 1993) to find the length of the 


vector: 


WsP= [it@Pae (2.8) 


Inner products of two functions f(z) and g(x) also can be produced in the same 


manner: 
(F(¢),9(2)) = 7 9° = [> F(a).o* (wae, 


where f (z) and g(x) are defined between 0 < x <_X. 


C. ORTHONORMAL BASES 


Any two vectors are called orthogonal if they are perpendicular to each other. 
A basis Vj, Vo,...V, is called orthonormal if all the elements of V;,V9,...,V, are or- 
thogonal to each other, and |/ V; |} = 1,7 = 1,2,...,n (Strang, 1976). 


V1, Vo,.--, Vn is said to form an orthonormal basis set if the inner product is 


0 when 1437 
1 when i=]. 


(Yas =| 


D. FOURIER TRANSFORM 


The Fourier Transform of a square integrable function f (t) is defined by 


Fw)= [fe *** =(F(, e), (2.4) 
which is called the Fourier Analysis Formula. The inverse Fourier Transform is given 
by | 

1 — put 
f(i)= — | F (w) etd, (2.5) 


which is the Synthesis Formula. 














E. SAMPLING 


The Sampling Theorem allows us to convert the analog signal to a discrete 


form for discrete signal processing. The main idea is to take samples from a signal at 


such a rate that the sequence of samples uniquely defines the original analog signal. 


Let fr (t) denote the signal obtained by weighting the periodic sequence of 


delta functions spaced T seconds apart by the sequence of numbers { fr (t)}; then, 


the sampled signal is 


fr(t)= )) f (nT) 6(t- nT). 


m= —-CoO 


F. DISCRETE FOURIER TRANSFORM 


The Discrete Fourier Transform of a square summable function 
f (n), n =0,1,2,...N, is defined by 


N-1 


Elk = f (nea kn 
n=0 
Its inverse is given by 
1 lea 27 
fm=— >) Fk)ir™ 
N k=0 


G. SHORT TIME FOURIER TRANSFORM (STFT) 


1. Introduction 


(2.6) 


(2.8) 


It can be seen from Equation (2.4) that the Fourier Transform of a signal 


does not reflect the time-evolution of the frequencies (Vetterli and Kovacevic, 1995). 


Because it does not show how the frequencies vary with time, it is not possible to 


extract the transient anomalies in the signal, which is spread over all frequencies. 


Gabor (1946) proposed a time-frequency localization technique, which is known 


as the Short Time Fourier Transform (STFT). In this transform, the signal first is 


divided into short consecutive segments, and then the Fourier Transform of each seg- 
ment is computed. T’he problem with this transform is the poor time resolution if the 


transient phenomena is shorter than the time window. 
2. The Continuous STFT 


In order to obtain the STFT of the signal, the standard Fourier Transform is 


modified by inserting a time window g(t) in Equation (2.4) as follows: 


Fst) = [f@ g(t — rT) e ” dt. (2.9) 


Equation (2.9) is called the Short Time Fourier Transform (STFT). When the 
window g(t) is Gaussian, the STFT is called the Gabor Transform. 


3. The Discrete STFT 


We can obtain the discrete STFT by sampling the continuous STFT, Equation 
(2.9), on a uniform grid (Akansu and Haddad, 1992) such that 


w= hwo, 7 = 2%. 


Thus, we obtain the discrete STFT 


F(k,n) = [- f(t) g(t — mr) eno dt (2.10) 


(f(t), 7°" g* (t — nt)). 


| 


This is the inner product of the signal f(t) and the modulated window g(t). 
The inner product yields large values when f (t) has sinusoids at or near the frequency 
kwo. The inner product is small at the resolution cell kwo, nT) when f (t) has sinusoids 
at locations other than kw. 

The inner product F' (k,n) is calculated for each window location n7. Then 
the window is shifted to the next location, (n + 1)79, and F'(k,n) is calculated for 
that resolution cell, kwo, (n+1)7). This procedure generates a two-parameter family 
F (k,n), which can be plotted on a time-frequency grid as shown in Figure 2.1 (Akansu 


and Haddad, 1992). 











Frequency 





Figure 2.1. Time-frequency grid for the discrete STFT 


The STFT, Equation (2.9) can also be interpreted as the convolution of g(t) 
2 g(—t) with the modulated signal e~7*f (t) such that 


F(w, t) = 9(t) *e7*f (8), (2.11) 


where g(t) is the time reversal of the window g(t) and “x” denotes the convolution 


operation. The convolution can be shown schematically as in Figure 2.2. 


f()$—eCOO g(t) F (w,t) 


e jut 
Figure 2.2. Convolution to obtain STFT 


If wo in Figure 2.2 is discretized, the modulated filter bank is obtained. If the 


output of Figure 2.2 is also sampled at t = nto, the corresponding discrete STFT is 


obtained, as shown in Figure 2.3. 











F (0wo, nT) SF (0, 2) 


F (Ow, t) 
3 g (Z) TO 
e J0wot A 


Wo FB lwo, To F 1% 
50) F (lw, t) (lw, 279) (1,%) 


T0 
~FJlwot 
f(t) “ —— F (20, t) \ F (2u0, 270) 2 F (2,1) 
Nee a 


€ 
3 5 Leet. Petomined Pe) 
G 
TQ 
e J kwot 


Figure 2.3. The discrete STFT as the modulated filter bank 


H. WAVELET TRANSFORM (WT) 


The Wavelet Transform (WT) is founded on a set of specific basis functions, 
which are called “little waves” or wavelets (Young, 1993). They include short du- 
ration/high frequency and long duration/low frequency functions. Each element in 
the wavelet set is constructed from the same function, which is called the “analyzing 


wavelet” or the “mother wavelet.” 
1. Continuous Wavelet Transform (CWT) 


The Wavelet Transform of a function f (¢), with respect to the mother wavelet 


w(t) is given by 


We = 2 [#@ ve (—) dt, ae Rt, BER — (2.12) 
a 1 t—b : 


= F(t), Te ¥ 





Ws 


a 


where superscript “x” denotes a complex conjugate, and a and 6 are scaling and 
esa 1 
translation parameters, respectively. The normalization term —= assures that the 


a 
energy of the scaled mother wavelet is the same as that of the mother wavelet. 
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There are three conditions for a function to be a mother wavelet. It must be 
oscillatory, it must decay to zero, and it must integrate to zero (Young, 1993). 

As can be seen from the Equation (2.12), the wavelet family incorporates 
shifted and scaled versions of the mother wavelet, w(t). If we define the wavelet 
family as { wap (t) an, where 

1 t—b 


Wap (t) = a uy oom (2.13) 


then this enables us to rewrite Equation (2.12) as 


Ws (a, b) 


[fe va(e) a (2.14) 
(f (t), Wab (t)) : 


The Continuous Wavelet Transform is invertible, provided that the mother 
wavelet 7 (t) satisfies the admissibility condition (Young, 1993). A wavelet is admis- 
sible if the admissibility constant, Cy, satisfies the following inequality 

= 
Cy= | EE ee. (2.15) 
0 Ww | 


where W (w) is the Fourier Transform of the wavelet ~ (t). 
If the admissibility condition is met for the wavelet function, the inverse CWT 


exists, and is given by 


(oS 7 [ 7 / ~ = Wj(a, b) Was (t) dadb. (2.16) 


2. Discrete Wavelet Transform (DWT) 


The Discrete Wavelet Transform (DWT) can be obtained by sampling a,b in 
Equation (2.13) at 


a=aj’?, b=kbgag’. (2.17) 


Then the wavelet family takes the form 


Wye (t) = ag!” (ast — hbo), (2.18) 


where j,k € Z. If this wavelet family is complete in L? (R) for some choice of w(t), 
a, and 0, then the set {5x} are called frames (Vetterli and Kovacevic, 1995). Now 
f (t) in L? (R) can be expressed as 


F(t) = Dd) do die Vie (4), (2.19) 


where the wavelet coefficient dj, is the inner product 


aay 


| 


(F(t), de (8) (2.20) 
a!” { f(t) ¥* (abt — bo). 


Since the sets are frames not bases, they do not satisfy the Parseval’s Theorem. 
Also the expansion, Equation (2.19) is not a unique use of frames. 

A family of wavelet functions, {~;,} is said to be a frame if there exists 
0<A< B< oo, such that, for all f (¢) in L?(R), 


AIFOPS ED I Wm FO) PS BUFO, (2.21) 


where A and B are frame bounds (Vetterli and Kovacevic, 1995). 

When the two frame bounds are equal, the frame is called a tight frame. If 
w;, 18 a tight frame, with a frame bound A = 1, and if || ,, ||?= 1 for all j,k, then 
jx forms an orthonormal basis (Vetterli and Kovacevic, 1995). 


Since #;, is a tight frame, Equation (2.21) becomes 


ms Il (vin, fC) P= AU FO IP. (2.22) 


The reconstruction formula can be derived from the Equations (2.19), (2.20), and 
(2.22) as 


f@= rem Da ies F( t)) Dyn. (2,98) 


10 











When A = 1, Equation (2.23) is equivalent to Equation (2.19). 

Sampling a and 0 in CWT, Equation (2.13), as stated in Equations (2.17) and 
(2.18), makes the translation steps, b = kbgap J ; proportional to the time scaling, 
a = aj’. Accordingly, large scale values correspond to good time resolution, whereas 
small scale values correspond to poor time resolution. 

The uncertainty principle states that, for a given transform pair, f (t) @ F (w), 
the inequality 


l 
O10. > - (2.24) 


bo 


holds, where o; and o,, are defined as follows (Akansu and Haddad, 1992): 


pa tt lf@fae 2 Le lF we) Fw 
‘SF @)Pdt? "ff |Fw) Pde 


where the integration equation limits too are suppressed. 
As seen in Equation (2.24), there is a lower bound for the product o; o.,; thus, 
an enhancement in one domain results in degradation in the other domain. Figure 


2.4 shows the time-frequency tiling for the WT. 
' 


fo 
Figure 2.4. Time-frequency tiling for the WT 


3. Multiresolution Analysis 


a. Multiresolution Spaces 


A multiresolution analysis consists of a sequence of closed subspaces 


ll 





wV CG Vag VeC Vr e Vie (2:25) 


such that it has the following properties (Vetterli and Kovacevic, 1995): 


(a) Completeness: 


NV; = {0}, UV; =L?(R) jEZ (2.26) 
(b) Scaling Property: 


fe V; 0 Ff (2t) €E Vi41 (2220) 
foranyfunction f € L*(R) 


(c) Existence of Basis: 


There exists a scaling function ¢@ € V;, such that 


{bn (t)} = 27? o(Mt-—k) VjEZ (2.28) 


is an orthonormal basis for V5. 

Equation (2.25) is illustrated via a Venn diagram in Figure 2.5. The 
basic idea in multiresolution analysis is successive approximations. As can be seen 
from Figure 2.5, a signal in L* (R) can be approximated by a coarse representation 
V,;, which is spanned by $;,. The difference between the original signal and its 
approximation is the detail. This allows the vector space V; to be divided into two 


parts such as 


V; = Vj-1 8 W5-1 (2.29) 
Vet WW a 


where W;_, is defined as the orthogonal complement of V;_; in V;, and it provides 
the detail information to go from subspace V;_, to V;. The symbol “@” stands for 
direct sum, and it means that each element of V; can be written as the sum of the 


elements of W,; 1 and the elements of V;_1. 


i 











Figure 2.5. Equation 2.25 illustrated via a Venn diagram 
L? (R) also can be represented by the direct sum as 


L’(R)=V;,0 >> W;, (2.30) 
J 230 


as illustrated in Figure 2.6. 





Figure 2.6. Multiresolution representation of L? (R) 


b. Scaling and Wavelet Functions 

Equation (2.28) gives the scaling function, ¢;,, which spans the coarse 
approximation space V;, and Equation (2.18) gives the wavelet function, #,,, which 
spans the difference information space W,;. Together, they can represent any signal 
in V;41. It is also possible to interpret V; and W, as lowpass and highpass portions 


of the space Vj41. 
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Equations (2.18) and (2.28) imply that both ~;, and $;, are sets of 
translates of dilated functions of w; and ¢;, respectively. In Equation (2.28), 7 is 
a scale factor. When 7 > 0,¢;, narrows, V; gets wider accordingly, which results 
in finer detail. If 7 < 0,4,, stretches, V; gets narrower, and ends up with coarser 
resolution. | 

Let ¢ C Vo be a scaling function. The translates of the dilated ¢(t) 
span Vo. Thus, the scaling function ¢(t) can be represented as a linear combination 


of translations of ¢(2t) as 


g(t) = S— ho (k) 6(2t — k), KE Z, (2.31) 


k 
where the fo (k)’s form a sequence of real or complex numbers called scaling function 


coefficients. 
Let w(t) be an element of subspace Wy. The subspace W, itself is 


a subspace of V; spanned by ¢(2t). Therefore, w(t) can be written as a linear 


combination of translates of ¢ (2t): 


w(t) = Thi (k) o(2t — k), RET, (2.32) 


where the h, (k)’s are wavelet coefficients. 

It should be noted that since ¢(t) and w (¢) span two orthogonal spaces, 
ho (kK) and hj; (k) also should be orthogonal. Hence, orthogonality of Vo and Wo can 
be achieved by requiring that 


(ito (k), ha (K)) = 0. (2.33) 


This allows the wavelet coefficients, h; (k), to be related to the scaling coefficients, 
ho (k), such as 


hy (k) = (-1)* ho (N - 1 — &), (2.34) 


assuming NV is the even number of coefficients. Then, the inner product in Equation 


(2.33) becomes 
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(ag (Ka (K)) = SS ho (Bhi () (2.38) 
reas 
Nai 
= ho (k)(—1)"ho (N — 1 — k) 
= 0. 
Therefore, any function f (t) in L?(R) can be written as 
f)= s aide (t) + = . dyn Dir (t) , (2.36) 


k=— © j=l k=- 0 
where a, and d,;, are the Discrete Wavelet Transform Coefficients. T hey are defined 


as 


i 
en, 
hy 


Qk 


(t), dn) = [f(b (at (2.37) 
FOV) = [fF (ten (at, 


dix 


where ¢; (t) and #;, (t) are real functions. 

Equation (2.36) states that any function in L? (R) can be written as a 
linear combination of the scaling function at a fixed scale plus a linear combination 
of wavelet functions at higher scales. 

The scaling coefficients ho (k) should be chosen carefully to ensure that 
the resulting scaling function has compact support. For that reason, they must satisfy 
some necessary conditions given below (Newland, 1992). The development of these 
conditions are provided in Appendix A. 


(a) The sum of the scaling coefficients should be equal to two. 


k F 
(b) If a solution to Equation (2.32) exists, f ¢(t) dt = 0, and 
J | o(t) |?dt = 1, and an integer number of translations of ¢ () are orthogonal to each 
other, as defined by 


(#10). (0) = f $9 (t ~ Rat = 5(k) =| 
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then 


S- ho (n) ho (n ~— 2k) = 26 (k). (2.39) 


nr 


(c) The sum of the squares of the coefficients is equal to two. 
Let k = 0; Equation (2.39) then becomes 


% Hho (n) P= 2. (2.40) 


(d) The individual sums of the odd and even terms of hg (k) are equal 


to one. 


S1 + So — Zz 

where NV is even and 
N/2—1 
k=0 

and 

N/2-1 

Sp= > hy (2k + 1). (2.42) 
k=0 


c. Frequency Domain 


The DFT of hg (k) is defined as 


ECan ee 3 ho (e732 (2.43) 


Recall that Equation (2.31) can be interpreted as a convolution between ho (k) and 
o(2t). In this case, ho (k) can be viewed as an FIR digital filter impulse response 
and, accordingly, H (w) is the filter’s frequency response. 

If zero is substituted in frequency w in Equation (2.43), the d.c. re- 


sponse of the FIR filter is obtained as: 
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H (0) = Sv ho (k). (2.44) 


The right hand side of Equation (2.44) should be equal to two as the result of the 
necessary condition stated in Equation (2.38). Thus, the necessary condition, Equa- 
tion (2.38), is equivalent to requiring the FIR filter’s frequency response at d.c. to be 
two. 

If the integer translates of ¢ are orthogonal, then the constraint defined 


in Equation (2.38) is true, if and only if 


|H(w)/?+[A(w + a)? =4. (2.45) 


If equation (2.44) is viewed as an FIR filter, hg (k) is called a quadrature mirror filter 
impulse response (QMF) (Burrus and Gopinath, 1993). 


d. Filter Banks and Discrete Wavelet Transform 


Multiresolution analysis can be implemented by using a technique called 
Multiresolution Pyramid Decomposition or Mallat’s algorithm (Vetterli and Kovace- 
vic, 1995). 

Let f be a real function in V;.;. Since {23 * (27*7¢ — k)} spans Vj41, 


f can be represented as 


j+1 ; 
t)= lagen, 2°? 9(2t — k) = Dagan, oyaryp (t), (2.46) 
k k 
where 
agtiye = (F(t), OG 41k (2) y= [Fe ft p(27** t — k)dt. (2.47) 


Next, since V;,; = V;® W;, f can be expressed as the sum of two functions, one 


lying in V; and the other in the orthogonal complement W;. Hence, 


a? jeje (t) + y dj Wyx (t (2.48) 
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where the DWT coefficients a;, and dj, are given by 


| 


(d;n, f (t)) (2.49) 
(Vjk, I (t)) 


Qjk 


a; 


As mentioned earlier, the first term in Equation (2.48) is the coarse 
approximation, which is the lowpass portion of f (t), whereas the second term contains 
the detail information, which is the highpass portion of f (t). 

From Equations (2.31) and (2.32), we have 


# (t) = Sho (k) 4 (26 — ), 
v(t) = ofa (k) (28 — &) 


o(2)t —n) = Siho(k) o(2(2t - n) - &) (2.50) 


Sho (k) o(2?t1t — Qn ~ k) 


Let m = 2n + kork = m — 2n, which leads to 


¢(2t — n)=S~ho(m — 2n)h(27t4t —m). (2.51) 


Now, to obtain the DWT coefficients a,, and dj,, Equation (2.50) is substituted into 


Equation (2.49). Assuming all the values are real, 


aie = (jx, f @) | (2.52) 
[FO 2? ot - h) 

[fe 29/7 S~ho(m — 2k) 6(27+1t — m) 

Tho (m — 2k) 2-48 [™ Fe) 2g (aE — m) 


—&O 


= a S- ho (m = 2k) Q(j+1)m - 
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dig = 27"? Shy (m — 2k) agjaiym- (2.53) 


Equations (2.52) and (2.53) are convolutions followed by subsampling. 
These operations are illustrated in Figure 2.7, where ho (kK) and h, (k) are given by 


ho (k) = ho (—k) (2.54) 


hy (kK) = hy (—&). 


(541) 





Figure 2.7. One stage of multiresolution signal decomposition. 


Therefore, the Discrete Wavelet Transform Coefficients at scale j are obtained 
by convolving the coefficients at scale 7 + 1 with ho(—k) and h; (—k),.and then 
decimating to produce the expansion coefficients at scale 7. Figure 2.8 shows the im- 
plementation of Equations (2.52) and (2.53) for three levels. The notation LP stands 
for lowpass FIR filter whose weights are ho (—k), and HP stands for the highpass FIR 
filter whose weights are h, (—k). 
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Figure 2.8. Three levels of multiresolution analysis 


I. MEDIAN FILTERS 


A median filter consists of a sliding window. The center value in the window is 
replaced by the median of the values within the window. The median of n observations 


Zi,% = 1, 2,...,n is given by 


Ly ,b= 2v+] 

med (x;) = Wn, n - (2.55) 
2 (Zw) + Lv41)) N= 2v 

where z;;) denotes the center value of n observations (Pitas and Venetsanopoulos, 


1990). 
A median filter of size n = 2u + 1 is defined by the following input-output 


relation: 
vi = med aa 2+, D7, eq tae) Le Z ’ 


where z,; is the input sequence and y; is the output sequence. An example of median 
filtering of size n = 3 is shown in Figure 2.9. The input takes on three values: 0, 1, 
2. As can be seen in Figure 2.9, the output sequence is 3-valued and does not contain 


isolated jumps atn =2 andn=8. 
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Figure 2.9. (a) Input z to median filter. (b) Output y from median filter. 


Median filtering is a nonlinear operation. The main advantage of this kind 
of nonlinear filtering is its computational simplicity and robustness. Furthermore, it 
has good edge properties. As will be Seen in Chapter III, data processing, wavelet 
analysis, and signal identification are used to separate the signals according to the 
transmitter from which they originate. The data processing uses median filtering and 

extracts the envelope. The wavelet analysis uses a Daubechies wavelet in the Mallat 
| algorithm. The signal identification uses a distance criteria based on selected wavelet 


scales and wavelet coefficients. 
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lil. WAVELET SIGNAL PROCESSING 


A. INTRODUCTION 


Mallat’s Pyramid Algorithm was used to obtain the Discrete Wavelet Trans- 
form (DWT). All programs related to preprocessing and wavelet processing, as well 
as to signal identification are written in Matlab, a product of Mathworks, Inc. (1992), 
and are given in Appendix D. The signals used in this analysis do not lend them- 
selves directly to wavelet analysis. Therefore, the processing phase is divided into 
two parts: a preprocessing phase where the signals are exposed to linear and non- 
linear operations, and an identification phase, which will be explained in subsequent 
sections. 

The signals are turn-off/turn-on transients from push-to-talk radios. Nine 
transmissions from four different transmitters were used in the identification part of 
this thesis. 


B. PREPROCESSING PHASE 


1. Turn-On Signals from Push-To-Talk Radios 


The signals used in the analysis were collected and recorded by the Naval 
security Group Activity, Charleston, SC. Nine turn-on/off samples of each of four 
transmitters were recorded. All the radios were Motorola models. Each radio is 


identified by its model name and number and is tabulated below: 


Radio: | Model: 
Transmitter 1 (Tr1) Maxtrac 
Transmitter 2 (Tr2) Saber 
Transmitter 3 (Tr3) HT440 
Transmitter 4 (Tr4) Saber 


Tr2 and Tr4 are different radios of the same model. 

The carrier frequency of the radios was 138.525 MHz. The signals were filtered 
at 1 Mhz Bandwidth (BW), and then digitized with a sampling frequency of 5 Mhz 
and a center frequency of 1.075 Mhz. 

The recordings from each transmitter are numbered from one to nine. Figure 


3.1 shows the first turn-on signal of the four transmitters. It can be noticed that 
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the recordings differ from each other in duration as well as in the type of waveform 
produced during the transition part. The turn-off response was also examined, but 
did not lead to any identification clues and was, therefore, not included in the results 
section. | | 

The information useful for the identification of the transmitter lies in the en- 
velope, which is the low frequency component of the signals. Wavelet Transforms are 
not useful in analyzing low frequency, but are well-suited for short duration phenom- 
ena. ‘Thus, we need to transform the data into a form suitable for wavelet analysis. A 
preprocessing scheme is applied to the signals to transform them into a usable form, 
and to improve the signal-to-noise ratio (SNR), as presented in the next section. It 
should also be noted that the SNR’s of the recordings are not known, and a denoising 


process would definitely enhance the identification performance. 
2. Preprocessing Scheme 


There are four steps in the preprocessing phase: taking the envelope, median 
filtering, differentiating, and a final median filtering. Prior to taking the envelopes, 
the d.c. terms were removed. 100-point median filters were used for filtering the 
envelope, whereas a 50-point median filter was used after the differentiation. The 
sizes of the median filters were determined experimentally. The main reason for 
using median filtering was to avoid broadening the signal, which usually is caused by 
convolutional operations. 

The preprocessed signals for the first recordings from the four transmitters 
are shown in Figure 3.2. All of the four final pulse-shaped signal waveforms seem 
to be appropriate for WT analysis. The preprocessing scheme is illustrated step-by- 
step in Appendix C for the first signal of each transmitter. The final preprocessed 


pulse-shaped signals for all transmitters are included in Appendix C. 
C. IDENTIFICATION PHASE 
1. Reduced Set Representation 


One of the main drawbacks of the wavelet transform is the shift variance, 
since the wavelet coefficients of a signal and a shifted replica can be totally different. 
Therefore, a Euclidean distance measure (introduced by Aware, Inc., 1992) is used as 


the technique for the classification of signals. 
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Figure 3.1 Signal samples from four different transmitters 


(a) Signal 1 from Transmitter 1 (b) Signal 1 from Transmitter 2 
(c) Signal 1 from Transmitter 3 (c) Signal 1 from Transmitter 4 


Mallat and Zhong (1989) demonstrated that the maxima extracted from the 
modulus of the wavelet coefficients can be used to reconstruct the input signal. That 
is, the maxima of the modulus of the wavelet coefficients contain approximately the 
same amount of information as the original signal. Consequently, signal analysis can 
be performed based on the wavelet extremals, which form a reduced signal represen- 
tation. 

Thus, wavelet coefficients at each scale have been replaced by their extremal 
values. Hence, the reduced set is only nonzero where the scale has an extremal, and 
is equal to the original value at these locations. Wavelet scale coefficients, which are 


not extremals, are set to zero. 
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Figure 3.2 Preprocessed signal samples from four different transmitters 


(a) Signal 1 from Transmitter 1 (b) Signal 1 from Transmitter 2 
(c) Signal 1 from Transmitter 3 (c) Signal 1 from Transmitter 4 


2. Ranking/Matching Algorithm 

The first step in computing a pairwise distance measure is to rank the peaks 
of the sets to be compared. That is, the peaks at each scale are ranked by their 
amplitudes. The second step is to match the ranked peaks to form pairs. The peak 
with the highest rank in one set is matched to the peak with the highest rank in 
the other set. The next in rank is matched to that of the other set, and so on. To 
illustrate the ranking and matching scheme, let us consider f (k) and g(k), as shown 
in Figure 3.3. We first rank the peaks by amplitude. The amplitudes of the sig- 
nals in ascending order are -2, -1, 1, 2, 3, 5 for f (kK), and -1, 0, 0.5, 1, 2.5, 4 for 
g(k). Then, we match the peaks by taking corresponding ranked peaks and forming 
pairs. The peak with the highest rank in f (kK) is 5; the peak with the highest rank in 
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gure 3.3. Signals ‘3 be matched: 
(a) uate ee signal f(k (sme ) maxima of signal g(k) 


Table 3.1 Matched pairs of signals f (k) and g(k) 


(columns show matched pairs, ranked in decreasing order from left to right) 






g (k) is 4. We match these peaks with highest ranks first. Then, we match the next in 
rank, which are 3 from f (k) and 2.5 from g(k), and so on. The ranked and matched 
pairs are tabulated in Table 3.1. The matching scheme does not require the number 
of elements in both sets to be the same. When a peak in one set does not have 
a corresponding peak in the other, we insert a zero into the set that has a smaller 


number of peaks, and form the pairs by matching the remaining peaks and the zeros. 
3. Distance Measure 


The third step is to compute a distance measure for the matched pairs. The 
distance assigned to the pair is the sum over the Euclidean distances in each scale. 
Thus, 

d (a,b?) = 3s (Wi, (ah — 83)? ]? (3.1) 


(k,l) are locations of the matched peaks 
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where a? and }’ are the wavelet maxima at scale j, and a} and 0) are the values a? , 
at temporal locations k and I, respectively, and Wis is the weighting factor at scale 
j tor the relative distances between the corresponding coordinates of matched peaks. 


The weighting factor, Wi, is defined as 
Wi | nj > nj i (3.2) 


where n?, and nj are the coordinates of a} and bj, respectively. 

The similarity between two sets is described in terms of the sum of the Eu- 
clidean distances of amplitudes weighted by the square root of the relative distances 
between the corresponding coordinates. Basically, a small value | nj, —n} | and | a? —07 | 
imply a high degree of similarity between a’ and DB’. 

As can be seen from Equation (3.1), this method provides no penalty if the 
peaks happen to be at the same location, or the amplitudes of the peaks are equal to 
each other. The penalty weight, Wi. is determined by the separation of the peaks. 
Large separation between the matched peaks corresponds to large penalty factors. 
The distance measure also depends on the amplitude difference of the matched peaks. 
The distance measure is directly proportional to the difference in amplitude. It should 
be noted that even identical scale outputs with offset in time (due to signal delay) 


will have a nonzero weighting factor. 
4. Modification to Weighting Factor 


We now introduce a lineup procedure and a modified weighting term to offset 
the shortcoming of the technique, described in Section III.C.3. This is the most 
robust method for obtaining reliable distance measures, and is used to obtain the 
experimental results. | 

There are four different sets of signals. Each set belongs to a different trans- 
mitter. The method defined above was used to determine the transmitting source of 
the signal. For this purpose, a template for each transmitter set is needed to mea- 
sure the similarity with any signal of interest. ‘The objective is to apply the distance 
measure algorithm to identify the transmitters. 

Matching peaks by ranking preserves the order of the peaks when the signals to 
be compared are very similar (i.e., identical), which then causes the weighting factor 


in Equation (3.2) to be a constant denoted by W. For example, if we shift f (k) in 
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Figure 3.3 four units to the right and compare it to f (k), the weighting factor in W 
in Equation (3.2) is 4 for all matched peaks. Since the amplitudes of the matched 
peaks are the same, the distance measure defined in Equation (3.1) is zero, in spite 
of the nonzero weighting factor W. This is not the case when we add noise to one or 
both of the signals, because it will cause amplitude difference between matched peaks. 
Even though the signals are similar, a nonzero distance measure will be computed. 
Furthermore, this distance measure will be dominated by the constant weighting 
factor W. We need to reduce the weighting factor in such cases by eliminating the 
offset between the signals. To eliminate this problem, the maximum peaks of the two 
sets are aligned by removing the time shift offset between the signals. | 
The first signal of the four sets was chosen as a template. To improve the ac- 
curacy of the algorithm, every signal was aligned with each template before applying 
the distance measure, Equation (3.1). Hence, the relative distance between the max- 
imum matched peaks of the signal and templates is zero. Accordingly, the weighting 
factor, Wi, will be zero for the maximum peaks for all templates. Applying zero 
weight for the distance with the corresponding (i.e., proper) template makes sense, 
but not for the other templates. Therefore, the weighting factor is modified to be 
win { Ml Ezy 33) 
Now, a small distance measure implies more accurately that the signal belongs to the 


set from which the template came. 
5. Application of Distance Measure Algorithm 


The fourth, eighth, and sixteenth order Daubechies wavelet functions were 
used to compute the WT of the preprocessed signals. The distance measure algorithm 
was applied to measure the similarity between the wavelet coefficients of the signal to 
be identified and the templates. From experimental considerations, the eighth order 
Daubechies (Dau-8) wavelet functions gave satisfactory results, whereas Dau-4 and 
Dau-16 did not. The Dau-8 wavelet function, which was used as the mother wavelet, 


is shown in Figure 3.4. 
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Figure 3.4 Eighth order Daubechies wavelet function 


a. Transmitter 1 


The distance values at Scale 11 between the first four signals 
from Transmitter 1 and the templates are tabulated and shown in Table 3.2. As 
expected, the distance value with Template 1 is significantly less than that with the 
other templates, since Template 1 is the template for the signals from Transmitter 1, 


and small distance values imply high similarity. 


Scale 11 | Template 1 Template 3 


Signal 1 | 0 44.3073 14.7737 107.8438 
































Signal 2 | 3.0527 47.2379 16.3406 103.8587 
Signal 3 | 2.6195 40.5894 13.797 99.7881 
Signal 4 | 1.0806 42.9072 15.4847 105.308 


Table 3.2 Distance measures at Scale 11 for Signals 1 — 4 of Transmitter 1 


As can be seen from Table 3.2, the distance values, d, between signals from 
Transmitter 1 and Template 1 are well separated from those of the other templates. 
Zero value in distance indicates perfect match, and occurs between Signal 1 and 
Template 1. This is expected since Signal 1 was chosen as a template for the set. 

A thresholding technique can be introduced to automate the identification 
procedure by defining a threshold level for each template. For example,.one can 
compare the maximum value of Set 1 with the minimum values of other templates. 
Determining the threshold levels is not addressed in this study; however, comparing 


the minimum and maximum values of the template outputs may be quite useful to 
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evaluate the performance of the identification scheme. The maximum and minimum 
values quoted in the remainder of Chapter 3 are obtained from Appendix B. At Scale 
11 with Template 1, the maximum distance value is 5.3971, whereas the nearest 
minimum distance value of the other templates is 12.52. This minimum occurs with 
Template 3, which indicates that Transmitter 1 and Transmitter 3 have a somewhat 
similar output (i.e., turn-on behavior). The ratio of the minimum Template 3 output 
to the maximum Template 1 output is 2.319, which shows how well the Template 1 
results are separated from the other template results. This separation ratio is 2.381 
and 1.66 for Scale 10 and Scale 9, respectively. Lower scales were not useful in terms 


of the similarity measurement, and were not included. 
b. ‘Transmitters 2, 3, and 4 


Similarly, the application of the distance measure to the other 
three sets resulted in the following separation ratios. They show how well the signais 
of a particular set are separated from the other sets under worst case conditions. 
They are 3.56, 2.76, and 1.88 for signals from Transmitter 2 at Scales 11, 10, and 
9, respectively. The following separation ratios were obtained at Scales 11, 10, and 
9: 4.64, 1.98, and 3.73 for signals from Transmitter 3, and 1.75, 1.4, 1.24 for signals 
from Transmitter 4. Signal versus template output plots for all signal sets are given in 
Appendix C. For example, one can obtain minima, maxima, and standard deviation 
of the distance measures very easily. It should be noted that, for Transmitters 1, 2, 


and 3, nine signals were used while, for Transmitter 4, only six signals were usable. 
c. Probability of Identification 


Signal-to-Noise Ratios (SNRs) for all signals were estimated. As 
seen in Figure (3.1), the signal waveforms can be split into three parts. The off region 
is where there is no signal; the transition region is where there is a build-up from off to 
on state; and the on region is where the signal is at steady state. We can assume that 
the off region consists of noise only, allowing the computation of the noise power. 
Noise and signal coexist in the on region. The signal power can be computed by 
finding the power in the on region and subtracting the noise power. Hence, we can 
compute the SNR’s for all sets. An average SNR value of 32.38 dB, 40.87 dB, 38.58 
dB, and 31 dB was computed for Set 1, Set 2, Set 3, and Set 4, respectively. 
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Figure 3.5 Probability of identification (P;) of the signals 


(a) P; of signals from Transmitter 1 (b) P; of signals from Transmitter 2 


(c) P; of signals from Transmitter 3 (c) P; of signals from Transmitter 4 


Gaussian white noise was added to all sets to decrease the SNR levels. We 
repeated the experiment eleven times for different SNR values, and obtained the 
probability of identification versus SNR plots for the signals. Results are shown in 
F igure 3.5. Note that, when the signal of interest is from Transmitter 1 and its SNR is 
higher than 16 dB, then its identification probability is high, and the distance-based 
algorithm is reliable. The minimum SNR values required for a reliable classification 


at Scale 11 are 30, 23, and 30 dB for signals from Transmitter 2, 3, and 4, respectively. 
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d. Attempts to Reduce the Processing Load 


Using a reduced data set (i.e., taking the local extremal in dis- 
tance computation) reduced the computational complexity. In addition, two ap- 
proaches were attempted to obtain acceptable performance while minimizing the 
computational effort. In the first approach, the signal was decimated by using a 
Chebychev -2 filter after median filtering of the envelope. This reduced the num- 
ber of points and, therefore, the number of computations. The distance algorithm 
using the decimated data did not work well enough to classify the signals and was 
not pursued further. In the second approach, instead of matching all the maxima of 
the signal of interest and the template, only the highest four peaks from each were 
matched, and a distance measure was computed. Even the original data did not 
permit successful classification; hence, noisy versions of the signal were not tried. 

Matching of the peaks of the signal and template includes ranking. Therefore, 
it is highly affected by noise. To minimize the noise effect in the matching of peaks, 
point-by-point matching was used instead of ranking, but the results were not satis- 
factory. In summary, distance measures that use all extremals and match them by 


ranking are the most effective scheme in classifying the four signal sets. 
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IV. CONCLUSION 


A. CLASSIFICATION OF PUSH-TO-TALK COMMUNICATION . 
SIGNALS | | | 


The main objective of this study was to determine a wavelet-based algorithm 
designed to extract features to identify push-to-talk transmitters. Robustness, com- 
pact signal representation capability, and low computational complexity are the main 
advantages of the wavelet analysis, which was used for feature extraction. 

Push-to-talk communication recordings provided for this research have a com- 
mon feature: They all include a transition from the off-to-on state, as well as the 
on-to-off state. The turn-on transition phase is unique for each transmitter, and can 
be used for classification purposes. The recordings differ from each other by the wave- 
form in the transition region and the duration. The feature for the classification of 
the signals is contained in the signal envelopes. The turn-on part of the envelope is 
a transient and, hence, a wideband signal. 

The push-to-talk recordings were not in appropriate forms for WT analysis. 
Thus, the recordings were preprocessed to become usable for WT processing. Differen- 
tiation of the envelope of the signals was used to transform the data into pulse-shaped 
transients. Median filtering was applied to both the envelope and the differential to 
improve the signal-to-noise ratio. Median filtering was preferred since it does not 
broaden the transient features of the signal. 

A distance algorithm was introduced in this work. It was based on a Euclidean 
distance measure between the wavelet coefficients of two data set on a given scale. 
Decisions about the origin of the signal was made according to the distance measures 
between the signals and the templates, where each template represented a different 
transmitter. A small distance value implied that the signal belonged to the same set 
as that particular template. In its current form, the classification assumes that any 
signal of interest belongs to one of the four sets. 

The distance algorithm was applied to four different signal sets. The first 
recording in any of the sets was designated to be the template. Instead of using all 
the wavelet coefficients, just the local extrema were used. Using only edge points by 


taking the local extrema reduced the computational complexity in the algorithm. 
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The distance measure between the signal and the templates is the sum of all 
Euclidean distances between matched peaks of the signal and the templates. It also 
includes a penalty factor due to the relative square root distance between the matched 
ranked pairs and their difference in amplitude. Matching signal peaks to the template 
is performed by ranking. Before matching the pairs, the maximum peak of the signal 
is aligned with the maximum of the templates. This reduced the penalty weight for 
like signals that are not aligned in time. 

The WT of the signals was computed using Daubechies-4 (Dau-4), Daubechies- 
8 (Dau-8), and Daubechies-16 (Dau-16) as mother wavelets. Experimentally, it was 
found that Dau-8 worked best with this data set. 

Decimating the signal envelope before taking the differential reduces the num- 
ber of points and, hence, decreases the number of computations. The distance algo- 
rithm was applied to the decimated sequences, but the results did not allow classif- 
cation of the signals. 

We also considered using only the four highest wavelet coefficients at a given 
scale in the distance computation to reduce the computational complexity. However, 
simulations showed that the information cited was not sufficient. The distance algo- 
rithm described in Section III, using all the extremas, matching the peaks by ranking, 
and permitting a penalty according to Equation (3.3) allowed robust identification of 


the signal sets. 


B. RECOMMENDATION FOR FUTURE STUDIES 


The distance algorithm introduced in Chapter II] gave promising results in 
classifying the four signal sets provided for this research. Three important issues have 
not been addressed in this research: the template selection, a threshold technique, 
and incorporation of information from other scales. In this work, templates were 
chosen arbitrarily from the sets, and the distance algorithm was applied only to the 
four signal sets. When the signals to be identified are from these sets, the algorithm 
was capable of classifying the signals. However, if the signals do not belong to these 
sets, the algorithm will compute a distance to each of the templates, which can lead 
to misclassification. A thresholding technique and robust template selection should 


be the subject of further study. 
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Also, no attempt was made to combine information from several scales. Iden- 
tification was obtained by using just one scale. Typically, Scale 11 worked best. If the 
distance information from Scales 9 and 10 could be used, a more robust identification 


performance should be realized. 
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APPENDIX A. DEVELOPMENT OF NECESSARY 
CONDITIONS 


NECESSARY CONDITION 1: 


The condition given by Equation (2.38) can be derived by integrating both 
sides of Equation (2.31) (Burrus and Gopinath, 1993). 


/ $ (t) dt / Sho (k)b(2t—k)dt (A.1) 
_ Yoho (k) f 4 (2t - kat 


Let m = 2t — k, so dt = dm. Substituting into Equation (A.1) gives 


|| 


[swat = 5 hal . #(m)dm, (A.2) 
foltdt 1 
[o(m)\dm — 9 2 ho (F 


The left side of Equation (A.2) is equal to one, and so 
1 
~\ Ao (k) = 
y) a o ( ) i 


or 


\ > ho (kK) = 2. (A.3) 


NECESSARY CONDITION 2 


g@;, forms an orthonormal set for L? (R). Therefore, the inner product of ¢(t) 


with integer translate of itself is 
(6 (t), dx (#)) = | 6 (Ob tb) = 5 (h). (A.4) 
Substituting Equation (2.31) into (A.4) gives | 
[ SCX ho (n)ho (mm) ¢ (2t — n)9 (2 (t — k) — m)dt = 6 (k), 
3 Thao (mo (mm ) | $(2t-n)o )$(2(t — k) — m)dt = 6(k). 
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Let 2i=7, so dt = $ and 
2, ds ho (n) ho (m )= | o(r—n)b (7 - 2k - m)dt = 6(k). 
Since the integral is equal to 6 (n — m — 2k) from (A.4), then 
5 Dhaln m)5 (n —m — 2k) = 6 (k). (A.5) 
Let n = m+ 2k; then Equation (A.5) becomes 


Xho (n)ho (n — 2k) = 26 (). (A.6) 


NECESSARY CONDITION 3 
Let k = 0; then Equation (A.6) becomes 
Sho (n)ho(n) = 2, (A.7) 
De lho(n))? = 2. 


NECESSARY CONDITION 4 
Reproducing Equation (A.6) 

S— ho (n)ho (n — 2k) = 26 (k). 
Summing over k of both sides gives 


k nr 
Breaking Equation (A.8) into even and odd ordered terms, we obtain 
as ps ho (2m) ho (2k + 2m) + Sho (2k + 2m + 1)Ao (2m + »| =2 (A.9) 
k ™m m 


Rearranging terms results in 


m m k 


>. bs ho (2k + 2m) ho (2m) + > » ho (2k +2m+ 0) ho (2m+1)=2. (A.10) 
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Substituting Equations (2.41) and (2.42) into Equation (A.10), gives 
Si >> Ao (2m) + So S~ ho (2m +1) =2. 
Applying the same equations again into Equation (A.11) dieat 
Sy ho, So 
Equations 2.41 and 2.42 state that 
Sy + So = 2. 
Solving for 5; and 52 in Equations (A.11) and (A.12), gives 


Sy — S> ho (2n) = l, 


Sy = > hg (22 +1) = 1. 


4] 





(A.11) 


(A.12) 


(A.13) 


(A.14) 


(A.15) 
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APPENDIX B. DISTANCE MEASURES 


A. SIGNALS FROM TRANSMITTER 1 ARE THE INPUTS 





43 


TEMPLATE 1 
SCALE 6 SCALE7 SCALE8 SCALE9 SCALE10 SCALE 11 
SIGNAL 1: 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
SIGNAL 2: 0.0154 0.0189 0.1321 0.1166 0.5972 3.0527 
SIGNAL 3: 0.0119 0.0357 0.0512 0.1026 0.6823 2.6195 
SIGNAL 4: 0.0119 0.0090 0.0518 0.1234 0.3941 1.0806 
SIGNAL 5: 0.0096 0.0596 0.0811 0.1272 1.1857 4.2015 
SIGNAL 6: 0.0187 0.0144 0.1260 0.1907 0.5007 2.1553 
SIGNAL 7: 0.0088 0.0156 0.0644 0.2489 0.5075 3.4189 
SIGNAL 8: 0.0184 0.0353 0.0335 0.1883 0.5085 2.7106 
SIGNAL 9: 0.0108 0.0349 0.1020 0.4137 1.0445 0.3966 
TEMPLATE 2 
SCALE6 SCALE7 SCALE8 SCALE9 SCALE10 SCALE11 
SIGNAL 1: 0.2835 0.2481 0.5950 2.9081 9.6966 44.3073 
SIGNAL 2: 0.3811 0.2786 0.6446 2.4972 8.9884 47.2375 
SIGNAL 3: 0.3828 0.3385 0.4982 2.3074 8.6315 40.5894 
SIGNAL 4: 0.2978 0.2723 0.5886 2.3877 8.8995 42.9072 
SIGNAL 5: 0.3043 0.3478 0.4959 2.3277 8.1022 39.6375 
SIGNAL 6: 0.3253 0.2469 0.8060 2.3600 9.1666 43.6621 
SIGNAL 7: 0.2867 0.2623 0.5220 2.0454 9.1751 39.6036 
SIGNAL 8: = 0.3772 0.3061 0.5358 2.1604 8.8586 40.6839 
SIGNAL 9: 0.3756 0.3014 0.5133 1.9071 8.6736 39.2300 
TEMPLATE 3 
SCALE6 SCALE7 SCALE8 SCALE9 SCALE10 SCALE11 
SIGNAL 1: 0.0293 0.0478 0.1806 0.7911 3.2748 14.7737 
SIGNAL 2: 0.0195 0.0516 0.2315 0.7833 3.4095 16.3406 
SIGNAL 3: 0.0292 0.0873 0.1219 0.7441 3.3865 13.7970 
SIGNAL 4: 0.0147 0.0461 0.1442 0.7797 3.1756 15.4847 
SIGNAL 5: 0.0247 0.1011 0.1338 0.7385 2.8243 13.1612 
SIGNAL 6: 0.0143 0.0448 0.1829 0.6879 3.2091 13.4287 
SIGNAL 7: = 0.0215 0.0571 0.1285 0.8024 3.2493 13.7654 
SIGNAL 8: 0.0454 0.0721 0.1630 0.7261 3.3426 14.3787 
SIGNAL 9: 0.0333 0.0683 0.1201 0.7172 3.0344 12.5220 
TEMPLATE 4 
| SCALE6 SCALE7 SCALE8 SCALE9 SCALE10 SCALE1I1 
SIGNAL 1: 0.3094 0.5771 1.2600 4.7861 23.4442 107.8438 
SIGNAL 2: =: 0.3395 0.5918 1.3712 4.9459 21.3511 103.8587 
SIGNAL 3: 0.3230 0.6083 1.2067 4.6558 21.4529 99.7881 
SIGNAL 4: 0.3675 0.5963 1.2665 4.8759 22.1737 105.3080 
SIGNAL 4: 0.3766 0.6255 1.1815 4.7152 20.7366 101.1619 
SIGNAL 6: 0.3035 0.5860 1.3969 4.9315 22.9453 110.3185 
SIGNAL 7: 0.3733 0.6054 1.1867 4.3094 22.0092 95.2878 
SIGNAL 8: —(0..2758 0.5979 1.2183 4.4715 21.8996 99.7182 
SIGNAL 9: 0.3427 0.6121 1.1942 4.3671 21.8617 101.6282 


B. 


TEMPLATE 1 


TEMPLATE 3 


SIGNAL 1: 0.3185 0.2155 0.5388 2.9937 10.7275 92.0090 
SIGNAL 2: 0.1653 0.2228 0.4821 2.8090 10.7830 46.1492 
SIGNAL 3: 0.3053 0.3015 0.6483 2.4348 10.0656 48.1584 
SIGNAL 4: 0.2162 0.3174 0.7455 2.4423 10.9963 48.9938 
SIGNAL 5: 0.3320 0.3800 0.6067 2.6860 11.4619 99.0565 
SIGNAL 6: 0.2723 0.2287 0.6431 2.0230 11.5436 00.4853 
SIGNAL 7: 0.2483 0.3425 0.7296 2.9808 10.9103 59.0225 
SIGNAL 8: 0.3857 0.2846 0.8665 2.0800 12.4298 54.1137 | 
SIGNAL 9: 0.4016 0.3610 0.6024 2.8528 11.8614 57.9092 


TEMPLATE 4 








SIGNALS FROM TRANSMITTER 2 ARE THE INPUTS 


SCALE 6 SCALE7 SCALE8 SCALE9 SCALE10 SCALE1I1 


SIGNAL 1: 0.2835 0.2481 0.5950 2.5081 9.6966 44.3073 
SIGNAL 2: 0.1008 0.2111 0.5449 2.4071 9.9080 42.9127 
SIGNAL 3: 0.2231 0.2929 0.6255 2.2173 9.5509 45.2351 
SIGNAL 4: 0.1587 0.3089 0.7071 2.4259 10.2745 45.9598 
SIGNAL 5: 0.2397 0.3638 0.6821 2.6701 10.6592 50.7797 
SIGNAL 6: 0.1881 0.2625 0.6220 2.6232 10.7310 49.5033 
SIGNAL 7: 0.1906 0.3823 0.6591 2.3317 10.0235 50.3606 
SIGNAL 8: = 0.3259 0.2937 0.7425 2.0765 11.2478 49.7514 
SIGNAL 9: 0.3105 0.3787 0.6678 2.0760 11.1999 92.8772 
TEMPLATE 2 : 
SCALE6 SCALE7 SCALE8 SCALE9 SCALE10 SCALE 1I1 
SIGNAL 1: —_ 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
SIGNAL 2: 0.1504 0.1834 0.2414 1.1648 2.9675 8.3316 
SIGNAL 3: 0.0926 0.1289 0.2083 0.3845 2.0463 11.7426 
SIGNAL 4: 0.1542 0.1445 0.2714 0.2385 2.9824 10.2347 
SIGNAL 5: 0.1100 0.1513 0.1427 0.2236 1.4235 12.0260 
SIGNAL 6: 0.0919 0.0599 0.0953 0.2324 3.0988 3.5685 
SIGNAL 7: 0.1120 0.0864 0.2697 1.0350 0.4110 6.5288 
SIGNAL 8: 0.0429 0.0602 0.3798 0.2762 3.4554 2.7844 
SIGNAL 9: 0.0585 0.1117 0.2090 1.1790 3.0197 8.5858 


SCALE 6 SCALE7 SCALE8 SCALE9 SCALE10 SCALE 11 


SCALE6 SCALE7 SCALE8 SCALE9 SCALE10 SCALE II 


SIGNAL 1: 0.0688 0.3144 0.7227 2.4151 10.7828 46.0266 
SIGNAL 2: 0.2709 0.4742 0.7804 2.4630 12.0574 53.2623 
SIGNAL 3: 0.1439 0.2997 0.7299 2.3018 12.1275 51.6991 
SIGNAL 4: 0.2207 0.3734 0.7884 2.6920 12.1243 49.1001 
SIGNAL 5: 0.1785 0.2800 0.6635 2.7456 11.6775 49.5180 
SIGNAL 6: 0.1251 0.3715 0.8402 2.9717 12.1831 47.9240 
SIGNAL 7: 0.1568 0.2853 0.8519 2.6483 11.8194 52.1264 
SIGNAL 8: 0.1382 0.2722 0.8705 2.6883 11.7220 47.9139 
SIGNAL 9: 0.1242 0.3436 0.6856 2.4425 12.4913 590.0139 
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C. SIGNALS FROM TRANSMITTER 3 ARE THE INPUTS 








TEMPLATE 4 


TEMPLATE 1 
SCALE6 SCALE7 SCALE8 SCALE9 SCALE10 SCALE 1i 
SIGNAL 1: 0.0293 0.0478 0.1806 0.7911 3.2748 14.7737 
SIGNAL 2: 0.0258 0.0486 0.1886 0.7386 3.4439 16.0413 
SIGNAL 3: 0.0258 0.0488 0.1760 0.7199 2.9349 14.0412 
SIGNAL 4: 0.0243 0.0380 0.2365 0.8262 3.1571 15.9305 
SIGNAL 5: 0.0351 0.0413 0.2032 0.7480 3.3249 15.7214 
SIGNAL 6: 0.0181 0.0522 0.1437 0.6904 3.3154 15.6880 
SIGNAL 7: 0.0204 0.0407 0.1965 0.6858 3.0757 . 14.9835 
SIGNAL 8: 0.0211 0.0487 0.1396 0.7145 3.1253 14.8495 
SIGNAL 9: 0.0132 0.0385 0.1556 0.6820 2.6568 13.9337 
TEMPLATE 2 
SCALE6 SCALE7 SCALE8 SCALE9 SCALE10 SCALE11 
SIGNAL 1: 0.3185 0.2155 0.5388 2.5937 10.7275 52.0090 
SIGNAL 2: 0.2770 0.2549 0.6212 2.6437 10.7470 58.6234 
SIGNAL 3: 0.3123 0.2733 0.5398 2.9106 8.9611 54.0690 
SIGNAL 4: 0.3022 0.2595 0.6690 2.7903 9.6474 53.2208 
SIGNAL 5: 0.2388 0.2737 0.5984 2.6087 10.4443 58.9315 
SIGNAL 6:. 0.2715 0.2588 0.5428 2.5480 9.9815 52.8461 
SIGNAL 7: 0.2536 0.2180 0.5879 2.3883 10.0475 54.4567 
SIGNAL 8: (0.2530 0.2886 0.5516 2.6977 10.0750 50.2555 
SIGNAL 9: 0.2634 0.2249 0.4667 2.5582 7.2510 43.6765 
TEMPLATE 3 
SCALE6 SCALE7 SCALE8 SCALE9 SCALE10 SCALE11 
SIGNAL 1: 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
SIGNAL 2: 0.0145 0.0182 0.1264 0.1059 0.2600 2.5576 
SIGNAL 3: ‘0.0140 0.0574 0.0679 0.1825 0.5453 1.9164 
SIGNAL 4: 0.0087 0.0243 0.1422 0.0815 0.2531 0.5718 
SIGNAL 5: 0.0202 0.0205 0.1187 0.0708 0.2221 1.8509 
SIGNAL 6: 0.0173 0.0460 0.0522 0.1326 0.5117 1.8301 
SIGNAL 7: 0.0080 0.0164 0.2044 0.1456 0.2368 0.6173 
SIGNAL 8: 0.0202 0.0689 0.0991 0.1083 0.3890 1.7340 
SIGNAL 9: 0.0165 0.0355 0.0718 0.1247 1.3389 2.9981 


SCALE6 SCALE7 SCALE8 SCALE9 SCALE10 SCALE 11 


SIGNAL 1: 0.3651 0.5208 1.1902 5.0321 23.6918 115.5311 
SIGNAL 2: 0.3238 0.5603 1.3472 9.0435 23.8257 122.6571 
SIGNAL 3: 0.3244 0.5725 1.2440 9.3095 22.0900 119.5437 
SIGNAL 4: 0.3540 0.5590 1.3640 0.2345 22.9700 116.3537 
SIGNAL 5: 0.2868 0.5780 1.3037 9.0934 23.5456 124.5405 
SIGNAL 6: 0.3084 0.5596 1.2133 4.8809 23.0578 117.4324 
SIGNAL 7: 0.2942 0.5274 1.3532 4.6602 23.2007 118.0311 
SIGNAL 8: 0.2827 0.5685 1.2326 9.0956 23.2144 113.1257 
SIGNAL 9: 0.2940 0.5262 1.0909 4.9186 19.2677 104.3332 
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D. SIGNALS FROM TRANSMITTER 4 ARE THE INPUTS 
TEMPLATE 1 


TEMPL 


TEMPL 


SIGNAL 1: 
SIGNAL 2: 
SIGNAL 3: 
SIGNAL 4: 
SIGNAL 5: 
SIGNAL 6: 
SIGNAL 7: 
SIGNAL 8: 
SIGNAL 9: 


ATE 2 


SIGNAL 1: 
SIGNAL 2: 
SIGNAL 3: 
SIGNAL 4: 
SIGNAL 5: 
SIGNAL 6: 
SIGNAL 7: 
SIGNAL 8: 
SIGNAL 9: 


ATE 3 


SIGNAL 1: 
SIGNAL 2: 
SIGNAL 3: 
SIGNAL 4: 
SIGNAL 5: 
SIGNAL 6: 
SIGNAL 7: 
SIGNAL 8: 
SIGNAL 9: 


TEMPLATE 4 


SIGNAL 1: 
SIGNAL 2: 
SIGNAL 3: 
SIGNAL 4: 
SIGNAL 5: 
SIGNAL 6: 
SIGNAL 7: 
SIGNAL 8: 
SIGNAL 9: 


SCALE6 SCALE7 SCALE8 SCALE9 SCALE10 SCALE 11 


0.3094 0.5771 1.2600 4.7861 23.4442 107.8438 
0.3736 0.9836 3.0273 9.2057 38.3416 164.0387 
0.3331 0.6218 1.3112 5.3833 24.3415 113.3645 
0.3133 0.5338 1.3146 9.2073 24.1267 112.2255 
0.4374 0.5831 1.1027 5.0784 23.9667 107.1038 
0.3432 0.5722 1.2066 9.4055 23.0626 110.6970 
0.2560 0.4075 0.8567 3.8320 18.1956 83.7372 
0.2761 0.4673 1.3469 4.5214 28.6522 133.7564 
0.2699 0.6246 1.7105 7.0262 30.5762 136.1710 
SCALE6 SCALE7 SCALE8 SCALE9 SCALE10 SCALE 11 
0.0688 0.3144 0.7227 2.4151 10.7828 46.0266 
0.2104 0.6746 1.6970 4.8096 18.7294 72.9213 
0.1649 0.2812 0.8120 2.6968 11.2046 00.4846 
0.1370 0.2327 0.8829 2.6756 11.2253 07.7131 
0.1847 0.3027 0.5789 2.4825 10.9737 45.7909 
0.1039 0.2557 0.6717 2.8612 10.1952 53.3557 
0.1637 0.1727 0.5083 2.1826 6.2022 38.8847 
0.1866 0.2580 0.8640 2.9625 12.0382 54.5502 
0.1472 0.3361 0.9322 4.4756 14.4394 90.8276 
SCALE 6 SCALE7 SCALE8 SCALE9 SCALE 10 SCALE 11 
0.3651 0.5208 1.1902 9.0321 23.6918 115.5311 
0.4188 0.8757 2.3272 8.7571 35.0247 158.7921 
0.2974 0.5532 1.0881 5.0421 24.5287 120.8973 
0.3393 0.4954 1.1246 5.4472 24.5602 122.9673 
0.4295 0.5159 1.0414 5.1931 24.3348 114.7444 
0.3967 0.5028 1.0955 5.6244 23.1653 119.4670 
0.2582 0.3225 0.7155 4.1114 18.9958 94.1898 
0.2952 0.3749 1.1948 5.2218 29.6827 132.6012 
0.2697 0.5467 1.3014 6.8545 27.5281 133.7723 
SCALE6 SCALE7 SCALE8 SCALE9 SCALE10 SCALE 11 
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0.1732 0.5275 1.2713 3.4912 9.3186 30.3546 
0.0544 0.0466 0.3896 0.4135. 2.4614 6.5196 
0.0610 0.0694 0.4249 0.3959 2.6097 14.8425 
0.0986 0.0655 0.1514 0.4633 2.2296 6.8575 
0.0528 0.0506 0.0990 0.5903 3.7932 11.8557 
0.0814 0.1619 0.4059 1.7605 4.4264 22.1434 
0.0935 0.3151 | 0.6468 2.8861 13.5405 47.6043 
0.0776 0.3200 0.9673 3.8654 13.8608 49.1905 
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E. 


TEMPLATE 1 


SIGNAL 1: 
SIGNAL 2: 
SIGNAL 3: 
SIGNAL 4: 
SIGNAL 5: 
SIGNAL.6: 
SIGNAL 7: 
SIGNAL 8: 
SIGNAL 9: 


TEMPLATE 2 


SIGNAL 1: 
SIGNAL 2: 
SIGNAL 3: 
SIGNAL 4: 
SIGNAL 5:. 
SIGNAL 6: 
SIGNAL 7: 
SIGNAL 8: 
SIGNAL 9: 


TEMPLATE 3 


_SIGNAL 1: 
SIGNAL 2: 
SIGNAL 3: 
SIGNAL 4: 
SIGNAL 5: 
SIGNAL 6: 
SIGNAL 7: 
SIGNAL 8: - 
SIGNAL 9: 


TEMPLATE 4 


SIGNAL 1: 
SIGNAL 2: 
SIGNAL 3: 
SIGNAL 4: 
SIGNAL 5: 
SIGNAL 6: 
SIGNAL 7: 
SIGNAL 8: 
SIGNAL 9: 


SCALE6 SCALE7 SCALE8 SCALE9 SCALE 10 


0.0375 
0.0468 
0.0431 
0.0504 
0.0421 
0.0508 
0.0291 
0.0377 
0.0337 


0.1043 
0.1085 
0.0800 
0.1064 
0.0706 
0.0743 
0.1163 
0.1302 
0.1012 


0.2759 
0.3193 
0.2972 
0.3262 
0.3415 
0.3502 
0.3047 
0.2676 
0.2747 


0.7683 
0.7442 
0.6296 
0.7624 
0.6067 
0.8021 
0.6125 
0.5524 
0.5609 


1.3094 
2.0640 
1.3008 
1.7088 
1.5736 
1.5439 
1.2903 
1.3807 
1.5242 


SIGNALS FROM TRANSMITTER 1 ARE THE INPUTS: 
10 DB SNR REDUCTION 


SCALE 11 
3.3502 
3.6190 
3.5047 
3.4416 
5.0176 © 
5.0661 
3.8559 
3.7544 
0.8662 


SCALE6 SCALE7 SCALE8 SCALE9 SCALE10 SCALE11 


0.3553 
0.2227 
0.3561 
0.2451 
0.2911 
0.2959 
0.2792 
0.3114 
0.3268 


SCALE6 SCALE7 SCALE8 SCALE9 


0.0359 
0.0389 
0.0430 
0.0373 
0.0346 
0.0390 
0.0231 
0.0353 
0.0315 


SCALE6 SCALE7 SCALE8 SCALE9 SCALE 10 


0.3328 
0.2562 
0.3450 
0.3263 
0.3633 
0.3042 
0.3338 
0.3610 
0.3608 


0.2728 
0.2168 
0.2650 
0.2573 
0.2100 
0.2452 
0.3019 
0.2584 
0.2492 


0.1303 
0.0850 
0.0984 
0.1119 
0.0814 
0.0843 
0.1120 
0.1325 
0.1139 


0.5612 
0.4936 
0.5801 
0.5714 
0.5636 
0.5860 
0.4590 
0.5003 
0.5617 


0.6642 
0.6260 
0.5283 
0.7072 
0.5792 
0.6188 
0.5697 
0.5708 
0.5431 


0.3973 
0.3844 
0.3287 
0.3942 
0.3848 
0.4338 
0.3312 
0.3284 
0.3366 


1.3761 
1.2892 
1.1919 
1.3486 
1.2499 
1.3593 
1.2041 
1.2423 
1.2187 
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2.8099 
2.7765 
2.4908 
2.9680 
2.2985 
2.8187 
2.0446 
2.0022 
2.2212 


1.0658 
1.0643 
1.0871 
1.0760 
1.1143 
0.9025 
1.1282 
1.0011 
1.0887 


4.9135 
4.8175 
4.5966 
4.7544 
4.3789 
9.1565 
4.6390 
4.3976 
4.3478 


10.4138 
11.2037 
9.4202 
10.7271 
9.2627 
10.3953 
8.9549 
9.9918 
9.0768 


47.5912 
48.0249 
43.7669 
46.0365 
42.8591 
49.7207 
42.6443 
43.5053 
42.7932 


SCALE 10 SCALE 11 


3.8188 
4.3243 
4.0226 
3.9323 
4.0585 
3.6642 
4.3698 
4.0881 
3.4706 


22.4741 
22.7738 
20.6796 
23.1659 
20.7572 
22.9280 
20.0266 
21.1563 
21.0854 


15.6111 
16.0790 
16.3021 
16.3871 
15.6347 
14.2155 
17.4460 
16.2695 
13.8132 


SCALE 11 
107.4480 
102.1320 
100.6919 
106.4092 
97.1743 
115.2724 
96.3665 
98.0142 
100.2342 


F. 


TEMPLATE 1 


SIGNAL 1: 
SIGNAL 2: 
SIGNAL 3: 
SIGNAL 4: 
SIGNAL 5: 
SIGNAL 6: 
SIGNAL 7: 
SIGNAL 8: 
SIGNAL 9: 


TEMPLATE 2 





SIGNAL 1: 
SIGNAL 2: 
SIGNAL 3: 
SIGNAL 4: 
SIGNAL 5: 
SIGNAL 6: 
SIGNAL 7: 
SIGNAL 8: 
SIGNAL 9: 


TEMPLATE 3 


SIGNAL 1: 
SIGNAL 2: 
SIGNAL 3: 
SIGNAL 4: 
SIGNAL 5: 
SIGNAL 6: 
SIGNAL 7: 
SIGNAL 8: 
SIGNAL 9: 


TEMPLATE 4 


SIGNAL 1: 
SIGNAL 2: 
SIGNAL 3: 
SIGNAL 4: 
SIGNAL 5: 
SIGNAL 6: 
SIGNAL 7: 
SIGNAL 8: 
SIGNAL 9: 


SIGNALS FROM TRANSMITTER 2 ARE INPUTS: 
10 DB SNR REDUCTION 


SCALE6 SCALE7 SCALE8 SCALE9 SCALE10 SCALE 11 


0.2998 
0.1149 
0.2483 
0.1683 
0.2497 
0.1699 
0.2107 
0.3312 
0.3012 


0.0353 
0.1267 
0.0796 
0.1472 
0.0873 
0.1078 
0.1093 
0.0543 
0.0705 


SCALE6 SCALE7 SCALE8 SCALE9 


0.3378 
0.1748 
0.3212 
0.2274 
0.3371 
0.2661 
0.2570 
0.3939 
0.3870 


SCALE6 SCALE7 SCALE8 SCALE9 


0.1007 
0.2222 
0.1201 
0.2068 
0.1445 
0.1370 
0.1369 
0.1439 
0.1384 


0.2975 
0.2315 
0.3285 
0.2955 
0.3499 
0.3207 
0.4434 
0.3403 
0.4205 


0.1029 
0.1354 
0.1694 
0.1665 
0.1667 
0.1217 
0.1754 
0.1232 
0.1886 


0.2623 
0.2254 
0.3164 
0.3074 
0.3617 
0.2662 
0.3376 
0.3071 
0.3634 


0.2607 
0.2770 
0.3323 
0.3380 
0.2756 
0.3112 
0.3029 
0.2392 
0.3070 


0.6936 
0.4739 
0.6071 
0.6568 
0.6493 
0.6972 
0.8268 
0.8018 
0.7371 


0.3632 
0.2513 
0.4690 
0.4090 
0.2662 
0.4735 
0.6404 
0.5491 
0.4846 


0.6906 
0.4875 
0.7229 
0.7264 
0.6446 
0.7978 
0.9129 
0.9439 
0.8302 


0.8255 
0.6216 
0.8504 
0.8050 
0.6615 
0.9905 
0.9008 
0.8706 
0.9407 
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2.5272 
2.1895 
2.1489 
2.4697 
2.6333 
2.3188 
2.5254 
2.6353 
2.3549 


SCALE6 SCALE7 SCALE8 SCALE9 


1.0898 
0.9447 
0.6193 
0.8105 
0.9486 
1.5780 
1.3555 
1.1825 
2.8504 


2.3864 
2.3379 
2.2456 
2.4591 
2.0699 
2.9988 
2.2908 
2.4221 
2.3635 


2.0868 
2.0512 
2.2053 
2.4836 
2.1278 
2.0032 
2.7150 
3.2170 
3.0078 


9.2397 
9.2329 
9.2771 
9.9201 
10.4630 
9.9021 
10.3927 
10.2630 
10.5009 


SCALE 10 
3.2442 
2.8265 
3.2228 
1.2229 
2.9736 
5.9286 
8.2270 
6.1173 
9.8706 


SCALE 10 
10.3773 
9.7415 
9.7959 
10.1376 
11.1484 
10.7795 
10.6121 
10.1326 
10.7260 


41.5986 
42.4683 
44.6032 
46.1952 
46.5821 
43.4716 
45.1681 
47.7992 
47.2822 


SCALE 11 
9.5561 
8.8378 
9.6595 

14.2345 

13.1501 

21.2568 

18.1313 
8.1226 

19.6026 


SCALE 11 
48.5725 
46.8498 
48.5703 
49.7111 
51.7127 
48.9604 
47.4262 
54.1949 
49.2730 


SCALE 10 SCALE 11 


9.9048 
10.6348 
10.3880 
11.3920 
10.2407 
11.8533 
11.7296 
12.0501 
12.4255 


44.4863 
48.9550 
45.7483 
47.6053 
51.0403 
53.8810 
53.8566 
46.7004 
52.4257 








G. SIGNALS FROM TRANSMITTER 3 ARE INPUTS: 
10 DB SNR REDUCTION 





TEMPLATE 4 


SIGNAL 1: 
SIGNAL 2: 
SIGNAL 3: 
SIGNAL 4: 
SIGNAL 5: 
SIGNAL 6: 
SIGNAL 7: 
SIGNAL 8: 
SIGNAL 9: 


0.0193 


TEMPLATE 1 
SCALE6 SCALE7 SCALE8& SCALE9 SCALE10 SCALE I1 
SIGNAL 1: 0.0478 0.1274 0.5144 1.1785 4.0000 18.1846 
SIGNAL 2: 0.0365 0.1234 0.4115 1.2813 4.1553 18.5676 
SIGNAL 3: 0.0684 0.1219 0.3577 1.1853 3.8279 16.0824 
SIGNAL 4: 0.0523 0.1170 0.4138 1.1140 4.3548 18.4488 
SIGNAL 5: 0.0348 0.1336 0.5074 1.2930 3.9845 17.9401 
SIGNAL 6: 0.0550 0.1275 0.4415 1.1653 3.9053 18.0756 
SIGNAL 7: 0.0478 0.1284 0.4557 1.1081 3.6589 17.1620 
SIGNAL 8: 0.0531 0.0784 0.3782 1.1528 4.0583 17.6531 
SIGNAL 9: 0.0320 0.1208 0.3447 1.0956 3.6144 16.3529 
TEMPLATE 2 | 
SCALE6 SCALE7 SCALE8 SCALE9 SCALE10 SCALE 11 
SIGNAL 1: 0.2735 0.2305 0.8375 3.2786 12.6610 59.1120 
SIGNAL 2: 0.2489 0.2555 0.7085 3.4069 13.2948 64.9324 
SIGNAL 3: 0.2403 0.2307 0.6814 3.2335 13.1578 56.7576 
SIGNAL 4: 0.2543 0.2354 0.6180 3.1848 12.5049 61.7671 
SIGNAL 5:, 0.2546 0.1906 0.8137 3.4720 12.2753 62.2135 
SIGNAL 6: 0.2131 0.1717 0.8292 3.7610 11.4564 99.3982 
SIGNAL 7: 0.2271 0.2378 0.7183 2.9783 12.5854 62.4342 
SIGNAL 8: —_0..2246 0.2336 0.7436 3.1300 12.7692 62.2046 
SIGNAL 9: 0.2461 0.2098 0.6158 2.1632 10.8813 49.2295 
TEMPLATE 3 
SCALE6 SCALE7 SCALE8 SCALE9 SCALE10 SCALE I1I1 
SIGNAL 1: 0.0203 0.0815 0.6030 0.9418 2.9238 6.4598 
SIGNAL 2: 0.0201 0.0950 0.4735 1.0146 3.0909 9.5653 
SIGNAL 3: 0.0465 0.0833 0.3634 0.9566 2.7794 6.4284 
SIGNAL 4: 0.0362 0.0954 0.4588 0.7535 2.6418 7.8841 
SIGNAL 5: 0.0253 0.0951 0.5927 1.0052 2.3614 7.6033 
SIGNAL 6: 0.0485 0.0933 0.5202 1.0480 1.7231 5.9440 
SIGNAL 7: 0.0328 0.0992 0.5659 0.7845 2.2809 8.0012 
SIGNAL 8: 0.0354 0.0650 0.4643 0.8565 2.4123 8.5680 
SIGNAL 9: - 0.0851 0.3710 0.7626 2.3411 5.4612 


SCALE6 SCALE7 SCALE8 SCALE9 SCALE10 SCALE 11 


0.3287 
0.1893 
0.2463 
0.3308 
0.2908 
0.2010 
0.2183 
0.1728 
0.2523 


0.5182 
0.5261 
0.5399 
0.5394 
0.4859 
0.4756 
0.5214 
0.4902 
0.5030 


1.5587 
1.4451 
1.3549 
1.3779 
1.6405 
1.5749 
1.4519 
1.4167 
1.2223 
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5.8201 
5.8022 
5.8309 
9.0470 
6.0450 
6.3159 
9.0075 
5.4070 
5.0020 


20.1909 
26.7385 
29.7889 
25.6568 
29.2253 
24.4335 
29.3703 
29.9855 
23.1458 


122.5166 


129.8573 


119.4983 
127.0563 
127.3273 
122.7423 
126.9125 
126.2491 
109.1846 








H. - SIGNALS FROM TRANSMITTER 4 ARE INPUTS: 


10 DB SNR REDUCTION 


TEMPLATE 1 


SIGNAL 1: 
SIGNAL 2: 
SIGNAL 3: 
SIGNAL 4: 
SIGNAL 5: 
SIGNAL 6: 


TEMPLATE 2 


SIGNAL 1: 
SIGNAL 2: 
' SIGNAL 3: 
SIGNAL 4: 
SIGNAL 5: 
SIGNAL 6: 


TEMPLATE 3 


SIGNAL 1: 
SIGNAL 2: 
SIGNAL 3: 
SIGNAL 4: 
SIGNAL 5: 
SIGNAL 6: 


TEMPLATE 4 


SIGNAL 1: 
SIGNAL 2: 
SIGNAL 3: 
SIGNAL 4: 
SIGNAL 5: 
SIGNAL 6: 


SCALE6 SCALE7 SCALE 8 


0.3270 
0.3248 
0.3220 
0.4828 
0.3500 
0.2645 


SCALE 6 SCALE 7 SCALE 8 


0.0850 
0.1570 
0.1514 
0.2097 
0.1014 
0.1545 


SCALE6 SCALE7 SCALE 8 


0.3780 
0.3014 
0.3627 
0.4693 
0.3955 
0.2624 


SCALE 6 SCALE7 SCALE 8 


0.0328 
0.0725 
0.0789 
0.1255 
0.0563 
0.0804 


0.5990 
0.5290 
0.6044 
0.4838 
0.5793 
0.4183 


0.3153 
0.2157 
0.2960 
0.2221 
0.3442 
0.2012 


0.5205 
0.4644 
0.5140 
0.4223 
0.5023 
0.3113 


0.0486 
0.1947 
0.0741 
0.2016 
0.0770 
0.2145 


1.6784 
1.7946 
1.5821 
1.1432 
1.2490 
1.0173 


1.0099 
0.8610 
0.8447 
0.7012 
0.7740 
0.7999 


1.5290 
1.5550 
1.4117 
1.0334 
1.1421] 
0.9258 


0.4642 
0.2994 
0.3058 
0.2156 
0.2064 
0.7304 
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SCALE 9 SCALE 10 


5.8691 
5.8942 
6.1968 
5.2652 
5.8876 
4.2269 


SCALE 9 


3.3984 
3.0255 
3.6034 
2.9600 
2.8104 
2.0320 


SCALE 9 


9.9429 
5.8505 
6.1926 
5.4432 
5.9523 
4.2049 


SCALE 9 


1.0811 
0.9492 
1.0105 
0.7971 
1.9411 
1.9989 


25.0385 
23.5419 
25.1165 
22.8460 
23.7742 
18.5720 


SCALE 10 
13.1845 
13.4155 
12.9318 
10.1786 
13.7426 
8.7308 


SCALE 10 
24.3696 
22.9977 
24.7024 
23.3718 
24.0663 
18.3186 


SCALE 10 
4.1788 
2.1063 
4.5781 
1.8608 
4.4024 
7.3088 





SCALE 11 
118.3022 
117.8467 
118.5933 
112.4434 
113.3160 
84.9252 


SCALE 11 
62.8281 
03.4214 
50.1148 
54.1355 
49.5709 
34.0015 


SCALE 11 
121.4879 
121.4152 
123.3364 
117.9020 
119.0966 
87.9338 


SCALE 11 
23.1301 
11.3148 
9.8803 
17.1363 
11.8552 
29.4086 








APPENDIX C 


PREPROCESSED SIGNALS AND TEMPLATE OUTPUTS 
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Figure C.1 Preprocessing scheme to obtain an appropriate waveform for Signal 1 
of Transmitter 1. From top to bottom: Signal 1 of Transmitter 1 (d.c. removed); 
Upper envelope of Signal 1; Envelope after median filtering of size 100; Differential 
of median filtering output; Differential after median filtering of size 50. 
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Figure C.2 Preprocessing scheme to obtain appropriate waveform for Signal 1 of 
Transmitter 2. From top to bottom: Signal 1 of Transmitter 2 (d.c. removed); Upper 
envelope of Signal 1; Envelope after median filtering of size 100; Differential of median 


filtering output; Differential after median filtering of size 50. 
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Figure C.3 Preprocessing scheme to obtain appropriate waveform for Signal 1 of 
Transmitter 3. From top to bottom: Signal 1 of Transmitter 3 (d.c. removed); Upper 
envelope of Signal 1; Envelope after median filtering of size 100; Differential of median 


filtering output; Differential after median filtering of size 50. 
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Figure C.4 Preprocessing scheme to obtain appropriate waveform for Signal 1 of 
Transmitter 4. From top to bottom: Signal 1 of Transmitter 4 (d.c. removed); Upper 
envelope of Signal 1; Envelope after median filtering of size 100; Differential of median 


filtering output; Differential after median filtering of size 50. 
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Figure C.5 Preprocessed signals from ‘Transmitter 1. 
Top row (left to right): Signal 1; Signal 2; Signal 3 
Center row: Signal 4; Signal 5; Signal 6 
Bottom row: Signal 7; Signal 8; Signal 9 
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Figure C.6 Preprocessed signals from Transmitter 2. 
Top row (left to right): Signal 1; Signal 2; Signal 3 
Center row: Signal 4; Signal 5; Signal 6 
Bottom row: Signal 7; Signal 8; Signal 9 
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Figure C.7 Preprocessed signals from Transmitter 3. 
Top row (left to right): Signal 1; Signal 2; Signal 3 
Center row: Signal 4; Signal 5; Signal 6 
Bottom row: Signal 7; Signal 8; Signal 9 
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Figure C.8 Preprocessed signals from Transmitter 4. 
Top row (left to right): Signal 1; Signal 2; Signal 3 
Center row: Signal 4; Signal 5; Signal 6 
Bottom row: Signal 7; Signal 8; Signal 9 
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Figure C.9 Distance measures at the output of templates when signals from Trans- 
mitter 1 are the inputs. Horizontal axis shows the number of the signals. Each column 
represents a template. First row is for Scale 11; second row is for Scale 10; third row 


is for Scale 9. 
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Figure C.10 Distance measures at the output of templates when signals from ‘Trans- 
mitter 2 are the inputs. Horizontal axis shows the number of the signals. Each column 


represents a template. First row is for Scale 11; second row is for Scale 10; third row 


is for Scale 9. 
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Figure C.11 Distance measures at the output of templates when signals from Trans- 
mitter 3 are the inputs. Horizontal axis shows the number of the signals. Each column 
represents a template. First row is for Scale 11; second row is for Scale 10; third row 


is for Scale 9. 
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Figure C.12 Distance measures when signals from Transmitter 4 are the inputs. 
Horizontal axis shows the number of the signals. Each column represents a template. 


First row is for Scale 11; second row is for Scale 10; third row is for Scale 9. 
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Figure C.13 Distance measures when signals from Transmitter 1 are inputs and 
their SNR values are 10 dB lower than the original data set. Horizontal axis shows 


the number of the signals. Each column represents a template. First row is for Scale 
11; second row is for Scale 10; third row is for Scale 9. 
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Figure C.14 Distance measures when signals from Transmitter 2 are inputs and 
their SNR values are 10 dB lower than the original data set. Horizontal axis shows 
the number of the signals. Each column represents a template. First row is for Scale 


11; second row is for Scale 10; third row is for Scale 9. 
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Figure C.15 Distance measures when signals from Transmitter 3 are inputs and 
their SNR values are 10 dB lower than the original data set. Horizontal axis shows 
the number of the signals. Each column represents a template. First row is for Scale 


11; second row is for Scale 10; third row is for Scale 9. 
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Figure C.16 Distance measures when signals from Transmitter 4 are inputs and 
their SNR values are 10 dB lower than the original data set. Horizontal axis shows 
the number of the signals. Each column represents a template. First row is for Scale 


11; second row is for Scale 10; third row is for Scale 9. 
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APPENDIX D 
MATLAB CODES 


wtsignal.m 


Purpose: 
To transform the push-to-talk recordings into pulse shape forms which are suitable 


for wavelet analysis. 


Synopsis: 
wtsignal 


Description: 

The program preprocesses the data as mentioned in Chapter 3. The signal vector to 
be preprocessed is entered as the input. 

Created by Yalcin Payal, November, 1995. 


x=input(’enter the signal to be preprocessed = ’) 
x=x-mean(x); 

y=envelope(x); 

my=mdsmooth(y,100); 

dmy=diff(my); 

mdmy=mdsmooth(dmy,50); 

end 
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distance.m 


Purpose: 
To compute distance measure between the signal and the four templates. Signals are 


4096 points long. 


Synopsis: 


distance 


Description: 

Wavelet transform matrix of the signal and the templates are the inputs to the pro- 
gram. Program is written to work with the wavelet transform matrices produced 
by the subroutine “mapdn.m” The subroutine “mapdn.m” produces the magnitude 
square of the wavelet coefficients at each scale, “distance.m” needs the amplitudes 
of the wavelet coefficients. Thus, a modification to “mapdn.m” is essential prior to 
run “distance.m.” mapdn.m and other functions called by mapdn.m can be found 
in (Newland, 1992) or in (Pitta, 1995). dl, d2, d3, d4 are the outputs of the -pro- 
gram and they are distance measures between the signal and Template 1, Template 
2, Template 3, Template 4, respectively at each scale. 

Created by Yalcin Payal, November, 1995. 


clear, clg 

sl=input(’ enter the WT of the signal to be identified = ’); 
kk1=input( ‘enter the WT of the template 1= ’); 
mmi1=input(’enter the WT of the template 2= ’); 
ttl=input(’enter the WT of the template 3= ’); 
vvl=input(’enter the WT of the template 4= ’); 

for n=1:11 

t=length(kk1(1,:))/(2"); 

col=n+2; 

a=kk1(col,1:t:length(kk1(1,:))); 
b=mm1(col,1:t:length(mm1(1,:))); 
c=tt1(col,1:t:length(tt1(1,:))): 
f=vv1(col,1:t:length(vv1(1:))); 
s=sl(col,1:t:length(s1(1,:))); 


a=localext(a); 
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b=localext(b); 

c=localext(c); 

f=localext(f); 

d=localext(s); 

[temp1,i]=sort(a); 
[temp2,j|=sort(b); 
[temp3,k]=sort(c); 
(temp4,k]=sort(£): 
[x,m]=sort(d); 
shiftl=m(length(x))-i(length(a)); 
shift2=m (length(x))-j(length(b)); 
shift3=m (length(x))-k(length(c)): 
shift4—m (length(x))-q(length(£)); 
wl=|abs(m-i-shift1)]; 
w2=(abs(m-j-shift2)]; 
w3=|abs(m-k-shift3)]; , 
w4=|abs(m-q-shift4)]; w1(find(wi==0))=ones(1,length(find(wl==0))); 
w2(find(w2==0))=ones(1,length(find(w2==0))); 
w3(find(w3==0) )=ones(1,length(find(w3==0))); 
w4(find(w4==0))=ones(1,length(find(w4==0))) 
d1(n)=sum(sqrt(w1.*(temp1-x)’)); 
d2(n)=sum(sqrt(w2.*(temp2-x)?)); 
d3(n)=sum(sqrt(w3.*(temp3-x)?)); 
d4(n)=sum(sqrt(w4.*(temp4-x)”)); 

end 


} 
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localext.m 


Purpose: 

To take the local extrema of the signal. 

Synopsis: 

localext. 

Description: 

LOCALEXT extracts the local extrema of a vector. LOCALEXT (X) returns a copy 
of vector X with all non extreme values set to 0. 

[Y,K] = LOCALEXT (X) will also return the number of deleted samples of X. 
function [ou,k]=localext(sg) 


(c) Copyright 1994, by Universidad de Vigo 
Author: Sergio. J. Garcia Galan, e-mail: Uvi_WaveQ@tsc.uvigo.es 


l=length(sg); 
sg=sg(:)’; 
ou=zeros(1,1); 
ax=abs(sg); 

MX = max(ax); 
MX=MX-++1; 
ou(1)=sg(1); 
ou(1)=sg(1); 


k=0; 
for 1=2:]-1 
if (sg(i) > sg(i-1)) & (se(i) > se(i+1)) 
ou(i)=sg(i); 
k=k-+1; 
end 
if (sg(i) < sg(i-l)) & (sg(i) < sg(it+1)) 
ou(i)=sg(i); 
=k-+1; 
end 
end 
k=l-k; 
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envelope.m 


Purpose: 
To take the envelope of the push-to-talk communication signals. 


Synopsis: 
envelope . 


Description: 


[Y,M]=ENVELOPE (x) returns the envelope , Y , and modulation index, M, of the 
AM signal in vector x. 


Created by Dennis W. Brown, 1993. 

Naval Postgraduate School, Monterey, CA. 

May be freely distributed: not for use in commercial applications. 
Part of the SPC Toolbox (Brown, 1995). 


function [y,m] = envelope(x) 
y=(]; 
if nargin = 1, 
error(’envelope: One argument required...’); 
end 
if min(size(x)) = 1, 
error(’envelope: Input arg ”x” must be a 1xN or Nx1 vector.’); 
end; 
x= 5G); 
y = abs(hilbert(x)); 
mmax = max(y); 
mmin = min(y); 
m = (mmax - mmin) / (mmax + mmin); 
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mdsmooth.m 


Purpose: 
To filter the envelope of the push-to-talk communication signals and the differentia] 


of the envelopes by median filtering. 


Synopsis: 
mdsmooth. 


Description: 

[y] = MDSMOOTH(X,L) smooths the input vector X using amedian filter with a 
rectangular window of “L” samples. 

MDSMOOTH is implemented as a mex function on someinstallations. 


Created by LT. Dennis W. Brown 

Naval Postgraduate School, Monterey, CA 

May be freely distributed; not for use in commercial applications. 
Part of the SPC Toolbox (Brown, 1995). 


function [y] = mdsmooth(x,L) 


y=[; 

if nargin = 2, | 
error(’mdsmooth: Invalid number of input arguments...’); 

end; 


if min(size(x)) = 1, | 
error(’mdsmooth: Input arg ”x” must be a 1xN or Nxl vector.’); | 

end; 

<= x(-); 

Ns = length(x); 

y = zeros(Ns,1); 

x = |zeros(L/2-1,1); x ; zeros ; zeros(L/2,1)|; 

for k=1:Ns 
y(k,1) = median(x(k:k+L-1,1)); 

end 
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